Casimir Energy for a Coupled Fermion-Kink System and its stability 
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We compute the Casimir energy for a system consisting of a fermion and a pseudoscalar field 
in the form of a prescribed kink. This model is not exactly solvable and we use the phase shift 
method to compute the Casimir energy. We use the relaxation method to find the bound states 
and the Runge-Kutta-Fehlberg method to obtain the scattering wavefunctions of the fermion in the 
whole interval of x. The resulting phase shifts are consistent with the weak and strong forms of the 
Levinson theorem. Then, we compute and plot the Casimir energy as a function of the parameters 
of the pseudoscalar field, i.e. the slope of <f>{x) at x — (fi) and the value of (j>(x) at infinity (#o)- 
In the graph of the Casimir energy as a function of n there is a sharp maximum occurring when 
the fermion bound state energy crosses the line of E = 0. Furthermore, this graph shows that the 
Casimir energy goes to zero for y, — ^ 0, and also for /i — > oo when 9o is an integer multiple of ir. 
Moreover, the graph of the Casimir energy as a function of 8o shows that this energy is on the 
average an increasing function of Qo and has a cusp whenever there is a zero fermionic mode. We 
finally compute the total energy of a system consisting of a valence fermion in the ground state. 
Most importantly, we show that this energy (the sum of the Casimir energy and the energy of the 
fermion) is minimum when the background field has winding number one, independent of the details 
of the background profile. Throughout the paper we compare our results with those of a simple 
exactly solvable model, where a piece-wise linear profile approximates the kink. We find that the 
kink is an almost reflectionless barrier for the fermions, within the context of our model. 
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I. INTRODUCTION 



The Casimir effect arises from the distortion of the 
zero-point energy of a quantum field due to the pres- 
ence of non-trivial background fields or the imposition 
of non-trivial boundary conditions. This effect was first 
proposed by Casimir in 1948 [l|, |2j, when he predicted 
the existence of an attractive force between two neutral 
infinite parallel metallic plates in a vacuum, placed a few 
micrometers apart. Since Casimir's work, this effect has 
attracted much interest and many authors have calcu- 
lated the Casimir energy and the resulting force caused 
by the presence of non-trivial boundary conditions for 
various geometries such as parallel plates, cylinders, and 
spheres |ll[3T-il8j and other geometries |19l - [26| . Moreover, 
they have used many different regularization and analytic 
continuation schemes to remove the divergences. Some of 
these techniques are the heat-kernel method [13, HH , the 
Green function formalism [29l |. the mode number sum- 
mation method combined with the zeta function analytic 
continuation technique (30l - [33j , and the multiple scatter- 
ing expansions [34]. The first experimental attempt to 
observe this phenomenon was conducted by Marcus Spar- 
naay [35j in 1958. In this experiment two parallel metallic 
plates were used, and the results had a very poor accu- 
racy. In 1997, Steve K. Lamoreaux [3(1 H?) measured the 
Casimir energy with a high accuracy using a plate and 
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a metallic spherical shell, and this was the first success- 
ful experiment to verify the Casimir effect. Since then, 
many different experiments have been performed to mea- 
sure the Casimir energy for various geometries [38|-[44j]. 

As mentioned above, the zero-point energy can also be 
affected by the presence of non-trivial background fields. 
The background field is usually chosen to be a soliton. 
Also sometimes a very simple potential such as an elec- 
tric potential well is chosen as the background field. This 
simple choice renders the problems of vacuum polariza- 
tion and the Casimir energy exactly solvable [45]. The 
Casimir energy also contributes to the lowest order quan- 
tum correction to the mass of the soliton. Many authors 
use this correction for models containing solitons such as 
supersymmetric solitons [IB-EH . For most of the models 
with solitons, the problem is not exactly solvable. Some- 
times indirect methods such as the phase shift method 
which relates the derivative of the phase shift with re- 
spect to the momentum to the spectral deficiency in the 
continuum states are used to calculate the Casimir en- 
ergy [HHiiHil 

In this paper we calculate the Casimir energy for a sys- 
tem containing a Fermi field chirally coupled to a pseu- 
doscalar field which is prescribed and has the form of an 
isolated kink. In a previous work [62j we calculated the 
Casimir energy for a similar system where the soliton 
profile is approximated by a piece-wise linear function, 
and this renders the problem exactly solvable. The vac- 
uum polarization for that model has also been calculated 
[63| . For that problem we calculated the Casimir energy 
exactly and directly by subtracting the vacuum energy 
of the system in the absence from the presence of the 
disturbance which is the pseudoscalar field. Throughout 
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this paper we shall refer to that model as the simple ex- 
actly solvable model. However, the present model is not 
analytically solvable and we use the phase shift method 
to compute the Casimir energy. As usually happens, the 
presence of the disturbance, e.g. the kink, leads to the 
appearance of one or more discrete bound states and also 
changes the continuum wavefunctions as compared to the 
free case. These changes have many manifestations in- 
cluding induced vacuum polarization and Casimir energy 
of the system. We have previously investigated a simi- 
lar system where neither the Fermi field nor the pseu- 
doscalar field, with boundary values of a topologically 
non-trivial configuration, were prescribed. They were al- 
lowed to interact and the non-perturbative final results, 
i.e. the results beyond the first order "back reaction" 
ones, revealed that the actual solitary wave profile differs 
only very slightly from an isolated kink 64] . This prox- 
imity is one of our motivations to study the properties of 
the coupled fcrmion-kink system. 

One of the main purposes of this paper is to inves- 
tigate how the functional form of a pseudoscalar back- 
ground field (4>(x)) affects the properties of the coupled 
fermion-pseudoscalar field system. Some of these effects 
have already been investigated (see for example 66] ) . For 
example it is well known that changing the value of <f>(x) 
at spatial infinity (±6>o) affects the spectral deficiency in 
the continua, and this usually leads to what is called the 
adiabatic contribution to the induced vacuum polariza- 
tion [67|, [68[. Moreover, changing the spatial profile of 
the background field close to the center e.g. the value 
of the slope of 4>(x) at x = (^), changes the pattern 
of energy levels crossing E = 0. These crossings lead to 
what is usually called the non-adiabatic contribution to 
the vacuum polarization [(£3, 68]. The same changes in 
the spectrum of the system that lead to the induced vac- 
uum polarization, also affect its Casimir energy. In this 
paper we investigate the effects of the functional form 
of 4>(x) on the Casimir energy and stability of the sys- 
tem. In order to accomplish this we meticulously study 
and compare the properties of the system which has the 
kink as the background field with those of the simple 
exactly solvable model. For the comparison to be mean- 
ingful and refined, we choose both background fields to 
have the same 9q and \i. In other words, this comparison 
serves a dual purpose: we not only investigate the be- 
havior of the systems as a function of the parameters 9q 
and /i, but also investigate the difference between these 
two systems for the same parameters. The latter investi- 
gation allows us to explore the effect of the finer details 
of 4>{x) on the overall properties of the system. 

In section 2 we briefly explain how to find the contin- 
uum scattering wavefunctions of this system using the 
Runge-Kutta-Fehlberg method of order 6. Then, in sec- 
tion 3 we compute the phase shifts. Using the relation 
between the phase shift and the difference between the 
density of states in the presence and absence of the distur- 
bance, we can write an expression for the Casimir energy 
in terms of the phase shift. We then calculate and plot 



the Casimir energy as a function of the parameters of the 
pseudoscalar field, i.e. the slope of the pseudoscalar field 
4>(x) at x = (/i) and the value of the field at infinity 
(6*o )• In both cases we also show the results for the simple 
exactly solvable model, for comparison. In section 4 we 
add the Casimir energy to the energy of a system consist- 
ing of a valence fermion in the ground state, when there 
is a soliton as the background field and discuss the stabil- 
ity of the system. For each result obtained and displayed 
for our model, we also present the corresponding results 
of the simple exactly solvable model, for comparison. In 
section 5 we summarize and discuss our conclusions. 



II. THE SPECTRUM OF THE FERMION IN 
THE PRESENCE OF THE PRESCRIBED KINK 

We consider the coupling of a Fermi field and a pseu- 
doscalar field governed by the following Lagrangian 

C = tP (i^d^ - Me^ (lh5 ) V, (1) 

and we choose 4>(x) to be prescribed in the form of 
4>{x) — 77i/\/Atanh \mx/\/2\ which is an isolated kink. 
The parameters M and m denote the mass of the Fermi 
and pseudoscalar field, respectively. Our purpose is to 
calculate the Casimir energy of this system. For exactly 
solvable systems we usually compute the complete spec- 
trum of the fermion, including the bound states with 
their discrete energies and the continuum states, and 
then we can calculate the Casimir energy directly by sub- 
tracting the vacuum energy of the system in the presence 
and absence of the disturbance. However, the form cho- 
sen for the (j)(x) makes the Euler-Lagrange equation of 
ijj analytically unsolvable. Therefore, we have to use an 
appropriate numerical method to find the fermion spec- 
trum. In order to facilitate the numerical calculations, 
we take advantage of the solutions of the exactly solv- 
able model for choosing the initial values as well as com- 
parison purposes. The Lagrangian for the simple exactly 
solvable model mentioned earlier is also the one shown in 
Eq. ([1]). However, the form chosen for (j)(x) is a piece- wise 
linear which renders the problem exactly solvable. The 
4>(x) for this model is as follows 

{—6*o for x ^ — I, 
fix for — I ^ x ^ I, (2) 
+9 Q for I < x. 

This form for the pseudoscalar field along with the kink 
is shown in Fig. ((T|) . In this figure we indicate the param- 
eters /x, the slope of 4>(x) at x = 0, and ±0 , its values 
at the boundaries. 
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FIG. 1: The solid and dashed lines show <j>(x) for the kink 
and the simple exactly solvable model, respectively. The pa- 
rameters #o and fi are also shown in the figure. 



A. The bound state energies and wavefunctions 

The spectrum of a Dirac field coupled to a background 
field gets distorted as compared to the free case. These 
distortions can be observed as spectral deficiencies in the 
continua and also the bound states could appear due to 
the presence of the background field. 
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Choosing the representation 7 = oi, 7 1 = ia 3 and 
(T2 for the Dirac matrices, the Dirac equa- 
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tion for the Lagrangian of Eq. ([I} becomes 

iaidtip — &3d x ip — M [cos <f> (x, t) + ia 2 sin 4> {%, t)] ip = 0, 

(3) 



where tp 



i'2 



We define 
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iip 2 
iip 2 



(4) 



The equation obeyed by £(x, t) is 



tion for j/^s are as follows 

y 1 + M cos <j> (x) y 3 - Ey 2 - M sin (j) (x) y± = 0, (6) 

y 2 + M cos <f> (x) y 4 + Ey x + M sin (j) (x) y 3 = 0, (7) 

y' 3 + M cos (f> (x) y x + Ey 4 + M sin <p (x) y 2 = 0, (8) 

2/4 + M cos <f> (x) y 2 - Ey 3 - M sin c/> (x) yi = 0, (9) 

where prime denotes the derivative with respect to x. We 
have an additional (fifth) equation E' = 0. To solve this 
set of five coupled first-order ODEs we need five con- 
ditions on the initial and final boundary points of the 
domain of the spatial variable, which could be chosen 
to be just [0, 00) instead of (—00, +00), due to the in- 
variance of the Lagrangian under the parity. Then, we 
map the x-interval [0, 00) to [0, 1] by the transformation 
X = tanh(x). Hence, the two boundaries of X are X = 
and X = 1. For bound states we use the relaxation 
method and choose the conditions at these boundaries as 
follows. At X — we choose two conditions: one parity 
condition and one assigning a value to one of the ?/i(0)s. 
This value is allowed to change so as to normalized t/j( x )- 
At X — 1 we choose three conditions: three of the j/i(l)s 
are set to zero. From now on we rescale all the quantities 
of the problem with respect to the mass of the Fermi field 
(A/). 

Some examples of the bound states obtained from nu- 
merical results are depicted in Fig. @. The upper graph 
shows the bound energy levels of the fermion as a func- 
tion of \i at 9q = 7T, i.e. a soliton with winding number 1 
and the lower graph shows the bound energies as a func- 
tion of #0 for the slope /1 = 10. In both graphs the bound 
energy levels for the model with kink and the simple ex- 
actly solvable model are shown by solid and dashed lines, 
respectively. The ± signs refer to the parity of each of 
the bound states. Notice that the rather small differ- 
ence between the profiles of these two background fields 
produces considerable difference between the pattern of 
bound state energies. 



/ id x -E iMe^W (0\ f . 

In order to obtain the energies of the bound states of 
the fermion, we use a numerical method called the re- 
laxation method (In an earlier paper (frij we have con- 
sidered the bound states of this system in detail). This 
method is used for solving the boundary value problems. 
To solve TV real coupled first-order ODEs, we need N 
boundary conditions, some of them are to be imposed 
at one boundary and the rest at the other boundary. 
The relaxation method determines the solution by start- 
ing with a guess and improving it, iteratively. We sepa- 
rate the real and imaginary parts of the upper and lower 

components of f^^j^ as ^ 1 ( a; ) = Vi( x ) + W2(x) and 
£2(2^) = 2/3(2;) + iyi(x). Therefore, the equations of mo- 



B. The continuum scattering states 

Now we focus on the continuum states. We can 
obtain the eigenstates of the Hamiltonian of the sys- 
tem, which can be chosen to be also the parity eigen- 
states. These solutions have to satisfy the parity con- 
dition P£(x,t) = —a 2 t;(—x,t). However, we obtain the 
wavefunctions which describe the scattering of a plane 
wave incident on the scattering region from the left, in- 
stead. The reason for obtaining the scattering solutions 
instead of the parity eigenstates is as follows. In the 
scattering process the wavefunctions on the right-hand 
side, far enough from the scattering region, have the in- 
dependence in the form of e lkx . Thus, by extracting the 
factor e lkx from the solution, the remaining part of the 
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FIG. 2: The energies of the bound states of the fermion. The 
upper graph shows the bound energies as a function of /i at 
9o — n and the lower graph shows the bound energies as a 
function of 8q at /j, — 10. In both graphs the solid and dashed 
lines are for the model with kink and the exactly solvable 
model, respectively. The parity of each bound state is indi- 
cated on the graphs by ± signs. 

solution on the far right would approach a constant, inde- 
pendent of the spatial variable and this will simplify the 
numerical analysis considerably. We choose the wave- 
function for the scattering process to be in the form 

t ( ~.\ _ jkx ( m(x) + im(x)\ 

where in general ?7iS are real functions of x. Then, the 
equations of motion satisfied by rjiS would be in the fol- 



lowing form 

rj 1 + cos 4> (x) i] 3 - (E + k)rj 2 - sin <j> (x) 774 = 0, (10) 

rj 2 + cos 4> (x) 774 + (E + k)rji + sin <j> (x) 773 = 0, (11) 

r? 3 + cos 4> (x) X]x + (E — k)rji + sin <j> (x) ri 2 = 0, (12) 

i] 4 + cos 4> (x) i] 2 — {E — k)i] 3 — sin <j> (x) r\\ — 0. (13) 

This set of equations for the form chosen for 4>{x), 
which is kink, cannot be solved analytically. Thus, we use 
again an appropriate numerical method. Our purpose is 
to find the scattering solutions of the system, not the 
parity eigenstates. Therefore, we need to find the wave- 
functions for the whole interval (—00, +00). We solve 
this set as an initial value problem, using the so-called 
Runge-Kutta methods. In order to find the solutions 
with high accuracy, we use the Runge-Kutta-Fehlberg 
method of order 6. We take advantage of the simple ex- 
actly solvable model to determine the initial boundary 
values for solving the equations. We already have all the 
solutions of this model, including the wavefunctions for 
the scattering process [H, |65[. Since the form chosen 
for the pseudoscalar field is similar to kink as the spa- 
tial variable x tends to infinity, we can use the values of 
the scattering wavefunctions of this simple model as the 
initial boundary conditions for our model. Therefore, we 
start at x — +00 with the values of the scattering wave- 
functions for the exactly solvable model (after dropping 
the factor e lkx and up to the normalization factor) and 
go backward in the x-interval to find the values of the 
rjiS for all the mesh-points of the interval (—00, +00) by 
solving Eqs. (jlOtil^ji . Figure © shows rjiS for the case 
with the parameters #0 = 7r, /1 = 10, k = 3.0 and 
E = +y/k 2 + M 2 . This figure also shows r/iS of the sim- 
ple exactly solvable model, for comparison. Note that 
the oscillations are less pronounced on the left for the 
kink model. This indicates that there is less "reflection" 
from the kink, as is also evident from the graph of p(x) 
shown in Fig. It is worth mentioning that the kink 
is totally renectionless for the elementary bosons within 
the A(/> 4 theory. 



III. THE CALCULATION OF THE CASIMIR 
ENERGY USING THE PHASE SHIFT METHOD 



The Casimir energy for a system like ours, as is well- 
known, is given by the shift in the zero-point energies of 
fermionic modes due to the presence of the disturbance, 



and in general can be written in the following form 
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FIG. 3: The graphs of r/i(x)s as functions of the spatial variable x, for the parameters 8q = tt, p = 10, k — 3.0 and 
E = +\/k 2 + M 2 . The solid and dashed lines show the normalized rii(x)s for our model and the simple exactly solvable model, 
respectively. 
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FIG. 4: The graph of p(x) as a function of the spatial variable 
x, for the parameters 8o — tt, p = 10, k — 3.0 and E — 
+\/k 2 + M 2 . The solid and dashed lines show the p(x) for 
our model and the simple exactly solvable model, respectively. 
The amplitude of oscillations on the left is proportional to the 
probability of "reflection" and the amplitude on the right is 
proportional to the "transmission" probability. 



The first equality is the relation we derived in [62[ for 
the Casimir energy of a Fermi field in the presence of 
an arbitrary disturbance. The functions v 3 p (x) and Vk (x) 
are the normalized fermion wavefunctions for the nega- 
tive continuum states in the presence and absence of the 
disturbance, respectively. The functions Xi\>i {%) are the 
normalized fermion wavefunctions for the discrete bound 
states with negative energy and £^~ md denote the en- 
ergies of these negative bound states. In the last line of 



the above equation, the extra M/2 takes into account the 
contribution from the half-bound state of the fermion at 
E = -M in the free case. The factor (> sca (fc) - p s ca (k)) 
is the difference between the density of continuum states 
with negative energy in the presence and absence of the 
pseudoscalar field. 

In [62| we concluded that for the simple exactly solv- 
able model we can calculate the Casimir energy only from 
the negative states or only the positive states, or the av- 
erage of all of the states and the results are exactly the 
same in all cases. Since all the symmetries of the model 
with kink are the same as the simple exactly solvable 
model, the aforementioned argument is also true for the 
present model. 



A. The phase shift and Levinson theorem 

The difference between the density of the continuum 
states in the free and interacting cases can be written in 
terms of the scattering phase shift in the following form 

p(k)-p (k) = --^-5(k), (15) 

7T dfc 

where S(k) = 5 s k y (fc) + S sea (k), i.e. S(k) sums over the 
contributions from both positive and negative energies. 
This relation is also true for the sea and sky, separately. 
Therefore, the second term in the relation of the Casimir 
energy, Eq. (1141) , can be written in terms of the phase 
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shift, as follows 
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(16) 



In the second equality we have just subtracted a zero 
term from the original one. For the last equality we 
have integrated the expression by parts, since the final 
expression is more convenient for the numerical analy- 
sis. Therefore, we can compute the second term in the 
expression of the Casimir energy using the phase shift. 

Now, by comparing the coefficients of e lkx on the left- 
and right-hand sides of the scattering region, we can ob- 
tain the scattering matrix element, which is related to the 
phase shift as S(k) — e tS ^ . We know that for x — > —00 
the wavefunction would be a linear combination of e tkx 
and e~ lkx . Therefore, to obtain the coefficients of e tkx 
on the left-hand side, we define the wavefunction of the 
far left as follows 



ai + ia 2 \ e - ikx 
a 3 + tci4 



b\ + 162 
63 + 264 



Akx 



m + "72 \ e ikx 



(17) 

where t^s on the right-hand side are obtained from the 
numerical method for x — > — 00. Since this relation has 
been written for the left-hand side and far from the scat- 
tering region, all the elements et^ and bi are constants, 
independent of the spatial variable x. The wavefunc- 
tion in the left-hand side of Eq. (TT7)) has to satisfy the 
first-order equations of motion (Eq. ([5])) when x tends to 
minus infinity in this set of equations. Using the relation 
obtained from the substitution of this wavefunction into 
Eq. §5§ and using Eq. (JTTJ) , we can obtain all the constants 
ai and bi. It is important to note that in the numerical 
computations we can only approach x = —00 or equiv- 
alently X = — 1. This decreases the accuracy, when the 
slope of <fi(x) at x — is very small, since in this case 
the region of variation of the pseudoscalar field becomes 
extended. 

Now we check the consistency of the resulting phase 
shifts with the Levinson theorem. The weak form of this 
theorem for the Dirac equation is as follows [65[ 
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where N is the total number of bound states, includ- 
ing positive and negative ones, N t the total number of 



the threshold bound states at the given strength of the 
potential, and 7V t ° the number of bound states at zero 
strength of the potential, i.e. the free Dirac case. The 
strong form of the Levinson theorem relates the value of 
the phase shift at each boundary of the continua to the 
number of levels which have crossed those boundaries, in 
the process of building up the disturbance. This form of 
theorem can be expressed in the following form for k = 



5(0) = (N cxit - N cntcI ) it. 



(19) 



This relation holds for each of the continua, separately. 
For each continuum N cxit (N cntor ) is the number of the 
bound states that exit (enter) that continuum from that 
boundary (E = +1.0 or E = —1.0) as the strength of the 
potential is increased from zero to its final finite value. 
In this equation the threshold bound states, mentioned 
above, should be included in N CX i t and N entel , with the 
coefficient 1/2. Moreover, the strong form of the Levin- 
son theorem for k = 00 can be written in the following 
form 



5(oo) = (TVcntcr - Ncxit) 7T. 



(20) 



This relation also holds for each continuum, separately. 
For each continuum iV CX it (A^ntcr) is the number of 
the bound states that exit (enter) that continuum from 
k = 00 (E = +00 or E — —00) as the strength of the 
potential is increased from zero to its final finite value. 
As an example of our results, in Figs. ([5]) and ([6]) we plot 
the phase shift for our system as a function of k, for 
the parameters 9 = n and [i = 10. Figure ([5]) shows 
the <5 s k y (fc), i.e. the phase shift for the states with the 
positive energy +Vfc 2 + M 2 and Fig. (|6]) shows S sea ,(k), 
i.e. the phase shift for the states with the negative en- 
ergy — \/k 2 + M 2 . In both figures we also show the phase 
shift of the simple exactly solvable model with the same 
parameters, for comparison. It is easy to check that the 
phase shifts depicted in these figures are consistent with 
both the weak and strong forms of the Levinson theorem. 
In particular 5(E = ±00) = ±#0, which is consistent 
with the results of the adiabatic method of Goldstone 
and Wilczek l67l. 



B. The Casimir energy 

Now we can calculate the Casimir energy using Eq. (|T4|) 
and show the results in some figures. In Fig. J7]) we plot 
the Casimir energy as a function of (j, at #o = 7r, i.e. a 
soliton with winding number 1, for our model and for 
the simple exactly solvable model. As can be seen in 
both cases, there is a sharp maximum occurring when 
the bound energy level crosses the line of E — 0. This 
crossing happens for a larger value of \x for our model, 
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FIG. 5: The graphical representation of <5 s k y (fc), for fj, = 10 
and 9o — n. The results for the models with kink and the 
simple exactly solvable model are shown by the solid and 
dashed lines, respectively. For these parameters we can see 
from Fig. ([2| that two levels have exited from the sky and 
since iV° sky = 1 at E = +1.0, we expect (5 s k y (0)/7r = 3/2. 
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FIG. 6: The graphical representation of 5 aea ,(k), for fi = 10 
and 8o = 7r. The results for the models with kink and the sim- 
ple exactly solvable model are shown by the solid and dashed 
lines, respectively. For these parameters we can see from 
Fig. ([2]) that no level has entered the sea and since iV t ° s ca = 1 
at E = —1.0, we expect <5 SGa (0)/7r = —1/2. 

based on the bound energy levels shown in the left graph 
of Fig. ([2) • For the model with kink the bound energy 
level crosses the line of E = at fi ps 3.821, while for 
the simple exactly solvable model this crossing occurs at 
a lower value of the slope, i.e. [i « 2.957. Also, the 
value of the Casimir energy is lower in the case of kink. 
The largest difference between the graphs of these two 
models occurs around the maximum, as is shown in the 
zoomed box of this figure. However, when ji — > or 
fi — y oo, both graphs tend to the same values. In both 
models the Casimir energy reaches the expectable value 
of zero when the slope of the soliton at x = decreases to 
zero, i.e. when the vacuum energy approaches that of the 
trivial vacuum, despite the residual non-trivial boundary 
conditions. Also, both graphs have the same limit when 
the slope of the soliton tends to infinity. This limit is 
zero at 9q = mr, i.e. when we have a proper soliton with 
winding number n. However, for other values of 8 the 
Casimir energy is in general non-zero when fi — > oo in 
both models (see Eq. (3.9) in [62l]). 
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FIG. 7: The graphical representation of the Casimir energy 
as a function of n, the scale of variation of the soliton, at 
#o = i"- The solid and dashed lines show the results for the 
model with the kink and the simple exactly solvable model, 
respectively. In the box we focus on the small values of /i to 
show the details of the maximum and the differences between 
the results of two models. 



In Fig. ([5) we present the Casimir energy as a function 
of 6*o for /.i = 10, by the solid and dashed lines for the 
model with kink and the simple exactly solvable model, 
respectively. As can be seen, the Casimir energy is, on 
the average, an increasing function of 9q for both models 
and there are two mild cusps in each graph. Comparing 
these graphs with the right graph of Fig. @, we conclude 
that these cusps occur when the bound energy levels cross 
the line of E = 0. 



1. 




-0.5 - 

-1. ' ' ' 

0. 1 77 UL 2jt 

2 2 

Bo 

FIG. 8: The Casimir energy as a function of Oo, at /j, — 10. 
The solid and dashed lines show the results for the model with 
the kink and the simple exactly solvable model, respectively. 
The cusps occur when bound state energy levels cross E = 0. 
The bound energy levels cross E = at 9o ~ 0.601-7T and 
8o ~ 1.801-7T, for the model with kink, and at 0o sa 0.5767T 
and Oo w 1.596-7T, for the simple exactly solvable model, in the 
interval ^ 0o ^ 2-7T. 



IV. STABILITY OF THE SOLUTIONS 

Finally, we consider the effect of the Casimir energy 
on the total energy of a system consisting of a valence 



8 



fermion in the ground state. The total energy for such a 
system consists of the Casimir energy and the energy of 
the valence fermion. Note that the energy of the valence 
fermion should not be added when this energy is nega- 
tive, since it has already been taken into account in the 
Casimir energy. The total energy is shown in Fig. ([9]). 
The upper graph shows this energy as a function of the 
slope of the pseudoscalar field (fi) when 9q = tt and the 
lower graph shows this energy as a function of the value 
of the pseudoscalar field at infinity (0 O ) when fi = 10. In 
both graphs the solid and dashed lines represent the to- 
tal energy for the model with kink and the simple exactly 
solvable model, respectively. As can be seen in the lower 
graph, for both models there is a minimum occurring at 
9q ~ tt, which corresponds to a soliton with winding num- 
ber one. This means that not only this configuration is 
energetically favorable, but also it is stable against small 
fluctuations in the parameters of the background field 
when this field is a soliton with a proper winding num- 
ber, as expected. Notice that there is no minimum in the 
upper graph. 




0.0 1 1 1 ' 1 1 1 ' ' 1 1 1 ' ' 1 1 1 ' 1 1 ' M 

5 10 15 20 




0.2 - 



0.1 ' ' ' ' 8„ 
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FIG. 9: The upper graph shows the total energy (the sum 
of the energy of a valence fermion and the Casimir energy) 
as a function of fi when 6q — tt. The lower graph shows the 
total energy as a function of 0o for the slope fi — 10. In both 
graphs the solid and dashed lines represent the total energy 
for the model with kink and the simple exactly solvable model, 
respectively. Notice that the system attains its lowest energy 

at 0Q = 7T. 



V. CONCLUSION 

In this paper we compute the Casimir energy for a 
Fermi field chirally coupled to a pseudoscalar field cho- 



sen to be prescribed and in the form of the kink. The 
equations of motion for this system are not exactly solv- 
able. Using the Runge-Kutta-Fehlberg method of order 
6 and taking the scattering solutions of an exactly solv- 
able model containing a fermion and a pseudoscalar field 
in the simple form of a piece-wise linear function, as the 
initial condition, we find the scattering wavefunctions for 
the fermion in the system with the kink. Comparing 
the graphs of the scattering probability for the system 
with the kink and the one with the simple background 
field shows that the kink is almost reflectionless for the 
fermion, within the context of our model. Then, compar- 
ing the values of these wavefunctions at the boundaries 
of the x interval, i.e. x — > ±oo, we obtain the phase shifts 
for the fermion's wavefunctions and plot them for the spe- 
cific set of parameters of the model. We check the con- 
sistency of the resulting phase shifts with the weak and 
strong forms of the Levinson theorem and conclude that 
this theorem is completely valid for our results. Then, 
using the relation between the derivative of the phase 
shift and the difference between the density of states in 
the presence and absence of the kink, we calculate the 
Casimir energy of our system and depict this energy as 
a function of the parameters determining the kink, i.e. 
its scale of variations (//) and its value at infinity (8q). 
We show that in the graph of the Casimir energy as a 
function of fi there is a sharp maximum occurring when 
the bound energy level of the fermion crosses the line 
of E = 0. Since this crossing occurs at a larger value 
of fi for the model with kink as compared to the sim- 
ple exactly solvable model, the maximum of the Casimir 
energy for the model with kink occurs at a larger value 
of fi. Moreover, this graph shows the expectable results 
at the limits fi — > and fi — > oo. For 9q = tt which 
describes a kink with winding number one, both limits 
are zero. The Casimir energy is in general an increas- 
ing function of 0q and is always positive and has a cusp 
whenever there is a zero fermionic mode. The first few 
cusps are local maxima. Finally, we compute the total 
energy of a system consisting of a valence fermion in the 
ground state. This energy includes the Casimir energy 
of the system and the energy of the valence fermion. We 
conclude that there is no preferable fi for the system con- 
sisting of the background field with winding number one. 
However, considering the effect of the changing the 0q 
for fixed fi, the system is more stable for the soliton with 
winding number one. 
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